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The recent BES Collaboration data on — > PV, particularly the isospin vi- 



\ dating mode '0' — > and finding of a finite number for B(ip' — > N 



enable us now to deal more precisely about the challenges to theory concern- 
■ ing this extraordinary and remarkable so called p — ir Puzzle of J ftp and ip 1 

. Decays. In terms of existing data, and deploying the simplest phenomenol- 

ogy measurement of i>' — > 7r + 7r and whether a finite number for the K* + K 
Q_i! mode might require a significantly larger accumulation of data, remain inter- 

esting questions. 



The strong suppression of ip' — > pir and K*K decays is well known and 
has been the subject for considerable experimental and theoretical activity. 
No evidence is seen for tp' — > pix in the BES group's 3.8 million tp' event 
sample, in striking contrast to the case for the J/ip, where J ftp — > pir is 
the dominant decay. This conundrum is referred to as the "p7r" puzzle and 
remains one of the most intriguing mysteries in quarkonium physics. The 
BES group has the largest event sample of ip' decays thus affording it the 
opportunity to undertake a systematic investigation of ip' decays to all of the 
lowest lying vector- plus pseudoscalar-meson (VP) final states. Of particular 
interest is the identification of the isospin violating mode ip 1 — > and 
the finding of a finite value for the branching ratio B(ip' — * K*°K ) which is 
likely an isospin violating VP case. These new numbers impel us by creating a 
"crisis" to deal in a systematic manner both theoretically and experimentally 
to pin down the essential nature of this so called p — it Puzzle of J ftp and 
ip' decays. Indeed interesting further but pedestrian tests appear feasible for 
the BES sample within a phenomenological framework, e.g. measurement of 
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ip' — ► ir + ir~ , cur] , pr/' and prj, though to reach for ip' to the charged K* + K~ 
VP mode might require a significantly larger data sample such as may be 
possible with the construction of a new main drift chamber MDC III and a 
future tau-charm factory at IHEP. 

On the theoretical side, there has also been much recent activity %. Li, 
Bugg, and Zou @ argued that final-state interactions involving the rescatter- 
ing of a\p and 02 into prr could be important and might interfere destructively 
in the case of ip'. The possibility of a destructive interference in ip' , though it is 
fortuitous, cannot be ruled out in view of Suzuki's finding H of the large long 
distance final-state interaction in J/ip based on existing data. However, this 
interference model does appear to have more assumptions than predictions! 
Brodsky and Kar liner §1 suggested that the decays J/ip,ip' — > pir proceed 
through the intrinsic charm component of the p wavefunction. They argue 
that the cc pair in the | udcc > Fock state of the p + has a nodeless wavefunc- 
tion which gives it a larger overlap with J/tp than ip'. This model dramatically 
challenges the assumption that charmonium states necessarily decay via inter- 
mediate gluons to exclusive low mass hadrons, deploying the analogy with the 
ss(4>) case found at LEAR where there is evidence in pp annihilation at rest 
of large violation of the OZI rule through gluon intermediary. However this 
paper is not yet characterized by experimentally checkable numbers for the 
charmonium system. On the subject of nodes in wave-functions, Pinsky % 
had proposed several years back that there is a node in the radial wave func- 
tion for ip', but not for J/ip, and that this node makes ip' — > pit a hindered 
Ml transition like J ftp — > r yrj c . However the BES measurements S ip' — > 77/ 
leads to the ratio 

B(ip(2S) -> jrf) m 

= WN> - W) = {) 

whereas Pinsky o relates the process ip' — > 77/ to the hindered Ml transition 
ip' — » 77/ c and predicts Q JV / = 0.002, which is well below the BES measured 
value. Finally there is the model of Chen and Braaten @ in which the cc pair 
is in a color octet 3 Si state for J/ip. 

Hou and Soni o developed an earlier suggestion of Freund and Nambu % 
that the decay J/ip — > pn is enhanced by the mixing of the J/ ip with a glueball 

that decays to pir. Brodsky, Lepage, and Tuan o (henceforth referred to as 
BLT) emphasized that J/ip — > pir violates the helicity selection rules of per- 
turbative QCD, and argued that the data requires the glueball O to be fairl 
narrow and nearly degenerate with the J/ip. It was however long recognized 
that the BLT model had a 'fly in the ointment' in the known isospin violating 
decay J/ip — > We shall discuss below this mode in some detail. Present 
data from BES constrains the mass and width of the glueball to the range 

1 mo - mj/ip I < 80 MeV and 4 MeV < T a < 50 MeV @0. As stressed by 
Chen and Braaten @ this mass is perhaps 700 MeV (or more) lighter than 
the lightest J FC = 1 glueball observed in lattice simulation of QCD with- 
out dynamical quarks (the "quenched" approximation) %. However, lattice 
calculations for the heavier glueball prediction may be less reliable than those 
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for the low mass J PC = ++ gluonium state, since as acknowledged by Bali et 
al.' 8 l lattice studies on the vector glueball are scarce and inconclusive, mainly 
because of the difficulties in constructing the corresponding lattice operators. 
Indeed recent QCD sum rule work also find m(l ) = 3.1 GeV(!), and a 3 
GeV 1 glueball is quite reasonable according to Brodsky @. 

A major problem for the Omicron O gluonium explanation of the pn 
puzzle, as well as other serious contending models which have an underlying 
assumption that hadron helicity conservation HHC holds at both J/ip and 
ip' (included here are the explanations of Brodsky and Karliner @ as well as 
Chen and Braaten fl which differ in details about the importance of DD 
channel and end point form factor H), is the known relatively large decay 
rate for J/ip -> n°u} !, the PDG branching ratio is (4.2 ± 0.6) x 10~ 4 , 
approximately three times larger than the J/tp — > ir + ir~ rate. Both of these 
1=1 decays are at first sight presumed to be due to electromagnetic decay 
(via a highly virtual ^(Q 2 ),Q 2 ^> 0) or gg^/, where the helicity properties 
of their amplitudes are identical to the strong decay J/ip — > ggg ^> pit %, 
and thus should satisfy the requirements of PQCD helicity conservation @. 
But there is no sign of suppression due to hadron helicity conservation for 
J/ip — > ujtt°\ One possibility is that there are additional qqg 1 = 1 resonances 
in the 3 GeV mass range which contributes to the ujtt channel (indeed the 
Omicron could be the 1 = member of a hybrid qqg state at 3 GeV). But 
this would destroy the elegance of Hou's argument ® on why a gluonium 
1™ should be at 3 GeV and in any case is a contrived band aid solution. 
Moshe Kugler suggested that the large then known J/ip — > ujit mode 
could be attributed to p — uj mixing. However the p — uj mixing model 
cannot give the large rate of J/ip — > ujit through J/ip — > p°ir —* lutt , 
because experimentally the rate of J/ip — > ujtt is as large as 0.1 times the 
rate of J/ip — > p°7r° 1111, whereas the effect of p — to mixing is only of order 
10 -2 (note that, e.g. the branching ratio of uj — > 7r + 7r~ is about 2%). Thus 
Brodsky presciently suggested even in 1989 H, that in any event it will be 
very important to compare these branching ratios for 7r + 7r~, and uj°tt° at the 
ip' and off resonance. 

One expects the J/ip and ip' mesons to decay to hadrons via 3g or for wzr 
via 7 or gg^ . In either case the decay proceeds via | ip(0) | 2 , where ^(0) 
is the wave function at the origin in the non relativistic quark model for cc. 
Thus it is reasonable to expect on the basis of perturbative QCD that for any 
exclusive hadronic final state h (including wvr°), the branching fractions scale 
like the branching fractions in e + e~ , to wit 

S = 0.UB(J/ip -» h)/B(ip' ^h)^ S e+e = 1. (2) 

Hence following up on Brodsky's concern B, it was argued that an in- 
triguing experimental measurement, if hadron helicity conservation HHC 
Theorem H is not applicable at both J/ip — ip' mass range, is to mea- 
sure ip' — > tt°uj at the level given by the 14% rule of Eq.d2D, i.e, 0.14 



xB(J/ip U7T°) ~ 0.6 x 10~ 4 . The BES data for B(ip' 



UTT°) iS 
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(0.40 ±0.10 ±0.06) x 10 -4 and hence taking into account experimental errors, 
the unr° decay is consistent with the 14% rule (2). It would appear that only 
one final test is needed to shut out the relevance of HHC for the J/ip/ip 1 
mass region, and this is to measure the ip' — > tt + it~ rate, which could be 
about three times smaller than ip' — > u>7r°, if we believe in the analogy with 
J /ip — > 7r + 7r~. Note also the stringent bound on B(J/ip — > 4>7r°) as well as 
the upper limit for B(ip' -> 0vr°) of < 0.66 x 10~ 5 (90% C.L.) in BES data § 
for these companion isospin violating decays. The situation here can however 
be understood in terms of the analysis of Haber &; Perrier @ , namely the re- 
duced branching ratio B(J/ip — > PV) = B(J/ip — > PV)/py(pv is momentum 
of vector meson in the center of mass) satisfy 



B(J/jj -> 7T 

B{J/ip -» 7r°o;) 



1 - (2)5 tan6»y 



tanfly + (2) 



for nonet symmetry. If (f>— oj are assumed to be ideally mixed as well [tan 6y = 
(1/2) a], then B{J/ip — > 7r 0) vanishes. A similar situation can be anticipated 
for £(V>' -» 7r°0). 

We must not however rush to judgment on the demise of HHC for J/ip/tp'. 
First, the experimental 'support' for HHC in J/ip — > lutt by Baltrusaistis et 
al. O is a qualitative one based on its steeper drop in form factor from q 2 = 
to q 2 = m 'j/ 1 p when compared with analogous irir form factor measured in 
J /if) — > ir + ir~ . Pakvasa £3 pointed out that it is known from the Sutherland 
Theorem based on PCAC that the isospin violating decay rp — > 7r + 7r~7r° 
vanishes for an electromagnetic intermediary, but rather is due to — m u 
current quark mass difference effect. In the present context JN) — > 
might have a dominant decay amplitude ~ (md — m u ) / m c = o/tt I from this 
'strong' effect not covered by HHC, and the electromagnetic intermediary 
7, gg 7 governed by HHC to be (a/7r)x [helicity suppression]. Perhaps we 
are seeing the former mechanism at work in the decays J/ip,ip' — > wk°. Our 
understanding of m u — rrid effect in charmonium decays is however clearly 
inadequate still. Though like Brodsky and Karliner @ (in a different context) 
we have dispensed with the 7*, gg 7, 3g contributions in ip' and J/ift decays 
to coir , it remains a mystery why for instance isospin violating ip' — > pir 
decay has not been seen through this mechanism U. Indeed the limit on 
ip' — > ujrj < 0.26 x 10 -4 (90% C.L.) when compared with if/ — > u>tt° rate 
may pose another example of this difficulty. Hence the present note is to be 
regarded as a stimulus for further theoretical study on the needed physical 
idea to complete our understanding. Setting aside this obstacle, it is amusing 
that a consistent picture can be sketched for all presently known ip' — > PV 
decays in which HHC is valid at both J/ip,ip' mass scales for the strong 
3g intermediary decays, based on the BLT model % supplemented by recent 
work S. 

Seiden, Sadrozinski, and Haber @ presented a general phenomenologi- 
cal parametrization of amplitudes for J /ip — > P + V in their Table IV. We 
adapt their analysis for ip 1 — > P + V using their notation. Hence g = strong 



(3) 
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singly disconnected SOZI amplitude with 3 gluons exchanged, e is the "elec- 
tromagnetic" amplitude which may actually arise from the (m^ — m u ) effect 
discussed above but expected to be comparable to the usual electromagnetic 
strength, r the doubly disconnected DOZI suppression factor breaking nonet 
symmetry. The SU(3) violation is accounted for by a factor (1 — s) for every 
strange quark contributing to g, a factor (1 — s p ) for a strange pseudoscalar 
contributing to r, a factor (1 — s v ) for a strange vector contributing to r, and 
a factor (1 — s e ), if relevant, for a strange quark contributing to e. One crucial 
ingredient to the analysis is the quark content of the rj and rf. We shall argue 
below that the BES data are in fact consistent with the assumption that 
the dominant part of the rj and rf wavefunctions consists solely of uu, dd, and 
ss. It has been known for sometime @ , that the cc contribution to r\ and rf 
is miniscule. Following , we write 

T) = X ri | uu + dd > /2 1/2 + Y ri | ss > 
rf = X' v \uu + dd> /2 1/2 + \ ss > 

and take their approximate values for the X's and Y's extracted from the two 
photon width of the rj, rf, to wit 

X v = 0.8, = -0.6 
X v , = 0.6, Y v , = 0.8 

For J/ip — > P + V , s is small ~ 10 - 20% of g and r ~ 0.15; we shall assume 
the same for ip' — > P + V . Indeed r, s, s p , and s v are all small numbers © for 
J/ip — > P + V. It seems reasonable, given the preliminary nature of the BES 
data for tp' — > P + V, to make the same small number assumption here. 
Where needed, and for definiteness, we take 

Sg — Sy — Sp — s — 0.15 

in rough accord with expectations in J/ip — > P + V @ . 

The analysis proceeds as follows. First tp' — > pir has not been seen, so 
each of p + 7r~, p°vr , p~ir + with amplitude g + e, must satisfy 

g + e = 0,OT g = -e. 

This leads to the remarkable conclusion that ip' — > pir, via 3 gluon strong 
(3g) decay is suppressed to the usual electromagnetic transition level ~ a/n 
in amplitude strength. Hence Brodsky-Lepage HHC Theorem appears 
to work for ip' — > pit strong (3g) decay, and furthermore Eq.(^) gives us a 
concrete measure of the magnitude of such helicity suppression. If we ignore 
small s, s e contributions, then for ip' — > K*K each of K* + K~, K*~ K + has 
amplitude 

g(l - s) + e{\ + s e ) ^g + e^0 

which is of course consistent with existing ip' — > K* + K~ upper limit 0. For 
K*°K°, K*°K° decay of ip', the amplitude for each is 



This is consistent with finite BES number for ip' — > K*°K , B {ip' — > 
K*°K°) = (0.84 ± 0.24 ± 0.16) x 10~ 4 fl. The central value is a little higher 
than the branching ratio for ip' — > which has also amplitude 3e. But 
within errors the two numbers are consistent; besides Eq. (^) is only approx- 
imate | g | may still be slightly larger than | e | leading ultimately to a finite 
branching ratio for ip' — > pir. To put in some numbers and give ip' — > K* + K~ 
a target estimate to shoot for, let us take s = s e = 0.15 for SU(3) breaking 
effects and use the expressions given on left hand side of (9) with the same 
phase space for charged and neutral members of K*K. We have [using the 
central value for B(ip' -> K*°K°)]. 

B(ip' -> K* + K~) = 1.008 x 10~ 6 

Assuming g + e = 0, and ignoring the small s p contribution (with phase space 
factor Pur)/pt v ' = 1-182), we have from Seiden et al. El Table IV adapted to 
ip' and Eq. (5) using the central value for the branching ratio B{ip' — > urj') @ 
= (7.9 ± 3.6 ± 1.5) x 10~ 5 

B(ip' -» urj) = (0.097) x 10~ 4 • 

This is still consistent with the BES upper limit 11 for B(ip' — > curj) < 0.26 x 
10 -4 . We shall not make estimates for ip' — > cjyq' from the known experimental 
number for ip' — > (j>r), because as seen from Seiden et al. £l the theoretical 
expressions here has the full gamut of parameters (g, s, s e , s v ,r, e, s p , X, Y). It 
would be easy but not particularly illuminating to obtain a branching fraction 
for ij)' — > 4>r]' compatible with the current experimental limit by adjusting the 
many parameters available in this case. The 'electromagnetic' transitions 
ij>' — > p°r},p°r}', where amplitudes @ are respectively ZeX^ and 3eX v /, can be 
related to the known ip' — > lott rate with amplitude 3e using Eq. (5). Taking 
the phase space factors into account bpo^/P^Ti- = 0-935 and Ppty/p^o = 
0.793] and again using the central experimental value for ip' — ► lott we have 

B(ip' fPrj) = 0.239 x W~ A ,B(ip' -> p%') = 0.115 x 10~ 4 • 

These results are quite consistent with BES values that B(if)' — > p°rj) = 
(0.21 ± 0.11 ± 0.05) x 10~ 4 and B(ij)' -» p°r]') < 0.3 x 10~ 4 (90% C.L.). The 
prediction of Eq. ( JTo) ) may require a larger sample of ip' than the 3.8 million 
currently available. This may be feasible with the construction of MDC III 
and a future tau-charm factory. We have not discussed the ij)' —* 777, 777' modes 
of BES because even though of the VP type, they involve an external 7 
which is not covered by the hadron helicity conservation HHC theorem of 
Brodsky-Lepage central to our discussion. 

To summarize, amongst the models proposed to explain the strong sup- 
pression in the pir, and K*K channel with some claim to a quantitative basis, 
the J/ip — O mixing model @il appears to survive the latest results, except 
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for isospin violating PV modes. This model requires the O mass to be quite 
near the J/ip and to have a substantial decay width to two-body VP me- 
son states, as particularly emphasized by BLT @. Although searches for this 
state in an e + e~ — > pir scan near the J/ij) mass @ and in -0' — > tt + tt^ pir 
decays have been negative, they have not provided very severe restrictions 
on the properties of the O. A better O search strategy may be needed. 
However, we have already seen from an analysis of the preliminary BES 
ijj' data for PV, that within this model framework, significant insights into 
hadronic physics have been obtained, in particular the strength of the HHC 
suppression to electromagnetic level at ip' and a new departure in the 
understanding of isospin violating decays of J/ip^tp' — » u>ir°, which however 
demands further study at a fundamental level. 

Remarks (i) The ifi' — > VT final states have been measured @ and found 
to be suppressed by a factor of at least three compared to the expectations 
of Eq. (2). Since hadronic VT decays conserve HHC @, some other mech- 
anism must be responsible for this (possibly mild) suppression in the BLT 
model Si and that of Chen and Braaten 6. This pattern can be explained 
by also taking into account the orbital-angular-momentum selection rule for 
exclusive amplitudes in perturbative QCD and is under study by Chen 
and Braaten. The Brodsky and Kar liner model can account for compara- 
ble suppression of both VT and VP modes since the tensor mesons here could 
have an appreciable intrinsic charm content. There is however a general be- 
lief that intrinsic charm content of low lying V, T mesons is small c3. (ii) 
Though the models of BLT &! as described here and Chen and Braaten o 
have predictions for ip' — ► lott, in order of magnitude agreement with exist- 
ing data 11 , they are based on the phenomenological treatment of Ref. [18] 
rather than a fundamental understanding of this isospin violating mode in 
the context of HHC, which remains a shared concern with the Brodsky and 
Karliner model % also. Furthermore, the need for Chen and Braaten % to 
invoke cc color octet 3 Si contribution to J/ip could be a dubious assumption, 
since Kroll 8 has arguments on why, for some exclusive J/ip decays at least, 
the color octet contribution is an 0{v 2 /c 2 ) effect and hence negligible (< 12% 
or so). Mahiko Suzuki 11 in his very thorough paper demonstrated that the 
existing data on J/^) -> ujtt , pix, K*K assuming that these decays are me- 
diated by conventional 7* and 3g respectively, would lead to the uncomfort- 
able situation of a large relative phase nearly 90° between the three-gluon and 
the one-photon amplitude. Hence he is pessimistic about extracting mean- 
ingful information on B-meson decay (presumably at the future B-factory 
facilities) on the parameters of the fundamental interactions here. Indeed 
he has gone further §1 to show that the large phase problem (and possible 
resolution of the puzzle) may rest with measurement of e + e~ — > 7 — > ir + ir~ 
and e + e _ — > 7 — > K + K~ off the J/ip without recourse to theory. However 
construction of MDC III and a future tau-charm factory may be needed to 
obtain this decisive measurement. We agree with Suzuki that a high pri- 
ority of future e + e~ facilities, e.g. the 5 x 10 7 J/if; events expected from the 
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BES upgrade/MDC III, is the high precision remeasurement of the J/ip decay 
branchings, particularly for p-K,K*K, and lott . (iii) Though it is hoped that 
eventually lattice gauge work would provide accurate prediction of the width 
of gluonium states like the Omicron, we note that the limit currently 611 set 
on Omicron width T < 50 MeV is already a severe one. At round 3.1 GeV a 
3g glueball cannot be too narrow U. Hence, we urge renewed effort at BES 
upgrade/MDC III to continue the search for the Omicron vector glueball by 
a scan of the J/ip resonance. Remember that the HHC theorem 1121 was de- 
veloped with very minimal assumptions from perturbative QCD. If it should 
prove to be invalid at the J/tp and tp' mass range, the validity of Eq. (2) can 
also be questioned. 
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ergy under Grant DE-FG-03-94ER40833 at the University of Hawaii-Manoa. 
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